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The  objective  of  this  article  is  to  present  a  detailed  review  of  the  literature  that  deals  with  the  appli¬ 
cation  of  CFD  in  the  design  of  solar  air  heater.  Solar  air  heater  is  one  of  the  basic  equipment  through 
which  solar  energy  is  converted  into  thermal  energy.  CFD  is  a  simulation  tool  which  uses  powerful 
computer  and  applied  mathematics,  to  model  fluid  flow  situations  for  the  prediction  of  heat,  mass 
and  momentum  transfer  and  optimal  design  in  various  heat  transfer  and  fluid  flow  processes. 
The  quality  of  the  solutions  obtained  from  CFD  simulations  are  largely  within  the  acceptable  range 
proving  that  CFD  is  an  effective  tool  for  predicting  the  behavior  and  performance  of  a  solar  air  heater. 
One  of  the  great  challenges  in  the  design  of  a  solar  air  heater  using  CFD  approach  is  the  selection  of 
appropriate  turbulence  model.  The  decision  about  a  suitable  turbulence  model  chosen  in  a  CFD 
computation  is  not  easy.  In  this  article  a  CFD  investigation  is  also  carried  out  to  select  best  turbulence 
model  for  the  design  of  a  solar  air  heater.  A  modern  CFD  code  ANSYS  FLUENT  vl2.1  is  used  to  simulate 
fluid  flow  through  a  conventional  solar  air  heater.  A  two-dimensional  flow  is  assumed.  The  influences  of 
the  five  different  turbulence  models  on  the  quality  of  the  obtained  results  are  tested.  It  appears  from 
the  performed  calculations  that  the  Renormalization-group  k-E  model  yields  the  best  results  for  two- 
dimensional  flow  through  conventional  solar  air  heaters. 
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Nomenclature 

/ 

friction  factor 

Nu 

Nusselt  number 

Ac 

surface  area  of  absorber  plate  (m2) 

Pr 

Prandtl  number 

cP 

specific  heat  of  air  (J/kg/K) 

p/e 

relative  roughness  pitch 

D 

equivalent  or  hydraulic  diameter  of  duct  (m) 

R 

roughness  function 

e 

rib  height  (m) 

Re 

Reynolds  number 

h 

heat  transfer  coefficient  (W/m2/I<) 

W/H 

duct  aspect  ratio 

H 

depth  of  duct  (m) 

I 

turbulence  intensity/intensity  of  solar  radiation 

Greek  symbols 

(W/m2) 

k 

thermal  conductivity  of  air  (W/m/K) 

(ra)e 

effective  transmittance-absorptance  product 

L 

length  of  duct  (m) 

t/th 

thermal  efficiency 

U 

inlet  length  of  duct  (m) 

P 

dynamic  viscosity  (N/s/m2) 

Li 

test  length  of  duct  (m) 

Pt 

turbulent  viscosity  (N/s/m2) 

L3 

outlet  length  of  duct  (m) 

P 

density  of  air  (kg/m3) 

m 

mass  flow  rate  (kg/s) 

(X 

angle  of  attack  (degree) 

A  P 

pressure  drop  (Pa) 

£ 

dissipation  rate 

P 

pitch  (m) 

CD 

specific  dissipation  rate 

Qu 

useful  heat  flux  (W/m2) 

5 

transition  sub-layer  thickness  (m) 

On 

useful  heat  gain  (W) 

K 

turbulence  kinetic  energy  (m2/s2) 

Ql 

heat  loss  from  collector  (W) 

r 

molecular  thermal  diffusivity 

Of 

heat  loss  from  top  of  collector  (W) 

r  t 

turbulent  thermal  diffusivity 

To 

fluid  outlet  temperature  (K) 

Ti 

fluid  inlet  temperature  (K) 

Subscripts 

To 

ambient  temperature  (I<) 

Tpm 

mean  plate  temperature  (I<) 

a 

ambient 

Tam 

mean  air  temperature  (K) 

am 

air  mean 

Tw 

wall  temperature  (I<) 

/ 

fluid  (air) 

Tm 

bulk  mean  temperature  (K) 

uL 

overall  heat  loss  coefficient  (W/m2/I<) 

V 

velocity  of  air  in  the  duct  (m/s) 

w 

width  of  duct  (m) 

pm 

plate  mean 

s 

smooth 

Dimensionless  parameters 

t 

turbulent 

w 

wall 

e\D 

relative  roughness  height 

1.  Introduction 

The  twenty-first  century  is  forming  into  the  perfect  energy 
storm.  Rising  energy  prices,  diminishing  energy  availability  and 
security,  and  growing  environmental  concerns  are  quickly  chan¬ 
ging  the  global  energy  panorama.  Energy  and  water  are  the  keys 
to  modern  life  and  provide  the  basis  necessary  for  sustained 
economic  development.  Industrialized  societies  have  become 
increasingly  dependent  on  fossil  fuels  for  myriad  uses.  Modern 
conveniences,  mechanized  agriculture,  and  global  population 
growth  have  only  been  made  possible  through  the  exploitation 
of  inexpensive  fossil  fuels.  Securing  sustainable  and  future  energy 
supplies  will  be  the  greatest  challenge  faced  by  all  societies  in  this 
century.  Due  to  a  growing  world  population  and  increasing 
modernization,  global  energy  demand  is  projected  to  more  than 
double  during  the  first  half  of  the  twenty-first  century  and  to 
more  than  triple  by  the  end  of  the  century.  Presently,  the  world’s 
population  is  nearly  7  billion,  and  projections  are  for  a  global 
population  approaching  10  billion  by  midcentury.  Future  energy 
demands  can  only  be  met  by  introducing  an  increasing  percen¬ 
tage  of  alternative  fuels.  Incremental  improvements  in  existing 
energy  networks  will  be  inadequate  to  meet  this  growing  energy 
demand.  Due  to  dwindling  reserves  and  ever-growing  concerns 
over  the  impact  of  burning  carbon  fuels  on  global  climate  change, 
fossil  fuel  sources  cannot  be  exploited  as  in  the  past.  Finding 
sufficient  supplies  of  clean  and  sustainable  energy  for  the  future  is 


the  global  society’s  most  daunting  challenge  for  the  twenty-first 
century.  The  future  will  be  a  mix  of  energy  technologies  with 
renewable  sources  such  as  solar,  wind,  and  biomass  playing  an 
increasingly  important  role  in  the  new  global  energy  economy. 
The  key  question  is:  how  long  it  will  take  for  this  sustainable 
energy  changeover  to  occur?  And  how  much  environmental, 
political,  and  economic  damage  is  acceptable  in  the  meantime? 
If  the  twenty-first  century  sustainable  energy  challenge  is  not  met 
quickly,  many  less-developed  countries  will  suffer  major  famines 
and  social  instability  from  rising  energy  prices.  Ultimately,  the 
world’s  economic  order  is  at  stake.  Approximately  one-third  of 
the  world’s  population  lives  in  rural  regions  without  access  to  the 
electric  grid,  and  about  half  of  these  same  people  live  with¬ 
out  access  to  safe  and  clean  water.  Solar  energy  is  unique  in  that 
it  can  easily  provide  electricity  and  purified  water  for  these 
people  today  with  minimal  infrastructure  requirements  by 
using  local  energy  resources  that  promote  local  economic  develop¬ 
ment  [1-5]. 

Solar  air  heaters,  because  of  their  simple  in  design,  are  cheap 
and  most  widely  used  collection  devices  of  solar  energy.  It  is  one  of 
the  basic  equipment  through  which  solar  energy  is  converted  into 
thermal  energy.  The  main  applications  of  solar  air  heaters  are  space 
heating,  seasoning  of  timber,  curing  of  industrial  products  and 
these  can  also  be  effectively  used  for  curing/drying  of  concrete/clay 
building  components.  A  conventional  solar  air  heater  generally 
consists  of  an  absorber  plate,  a  rear  plate,  insulation  below  the  rear 
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plate,  transparent  cover  on  the  exposed  side,  and  the  air  flows 
between  the  absorbing  plate  and  rear  plate.  A  solar  air  heater  is 
simple  in  design  and  requires  little  maintenance.  However,  the 
value  of  the  heat  transfer  coefficient  between  the  absorber  plate 
and  air  is  low  and  this  results  in  a  lower  efficiency  [6,7], 

CFD  is  a  science  that  can  be  helpful  for  studying  fluid  flow, 
heat  transfer,  chemical  reactions  etc.  by  solving  mathematical 
equations  with  the  help  of  numerical  analysis.  It  is  equally  helpful 
in  designing  a  solar  air  heater  system  from  scratch  as  well  as  in 
troubleshooting/optimization  by  suggesting  design  modifications. 
CFD  employs  a  very  simple  principle  of  resolving  the  entire 
system  in  small  cells  or  grids  and  applying  governing  equations 
on  these  discrete  elements  to  find  numerical  solutions  regarding 
pressure  distribution,  temperature  gradients,  flow  parameters 
and  the  like  in  a  shorter  time  at  a  lower  cost  because  of  reduced 
required  experimental  work  [8-13], 

A  lot  of  experimental  studies  have  been  carried  out  to  evaluate 
performance  of  solar  air  heater  but  very  few  attempts  of  CFD  investi¬ 
gation  have  been  made  so  far  due  to  complexity  of  flow  pattern  and 
computational  limitations.  With  the  development  of  computer, 
hardware  and  numerical  methodology,  applications  of  CFD  are  being 
used  to  carry  out  critical  investigations  in  the  field  of  solar  air  heaters. 
One  of  the  great  challenges  in  the  design  of  a  solar  air  heater  using 
CFD  approach  is  the  selection  of  appropriate  turbulence  model.  The 
decision  about  an  appropriate  turbulence  model  chosen  in  a  CFD 
computation  is  not  easy.  In  this  article  a  CFD  investigation  is  also 
carried  out  to  select  best  turbulence  model  for  the  design  of  a  solar 
air  heater.  The  main  objectives  of  this  article  are: 

1.  To  present  a  detailed  review  of  the  literature  that  deals  with 
the  application  of  CFD  in  the  design  of  solar  air  heater. 

2.  To  examine  the  predictive  ability  of  five  different  turbulence 
models  such  as:  Standard  k-e  model,  Renormalization-group 
fc-e  model,  Realizable  k-e  model,  Standard  k-co  and  Shear 


Stress  Transport  k-co  model  using  commercial  CFD  software 
ANSYS  FLUENT  v!2.1. 


2.  Thermal  performance  of  solar  air  heater 

Performance  of  any  system  represents  the  degree  of  utilization 
of  input  to  the  system.  It  is  required  to  analyze  thermal  and 
hydraulic  performance  of  a  solar  air  heater  for  making  an  efficient 
design  of  such  type  of  a  system.  Thermal  performance  concerns 
with  heat  transfer  process  within  the  collector  and  hydraulic 
performance  concerns  with  pressure  drop  in  the  duct.  A  conven¬ 
tional  solar  air  heater  (Fig.  1)  is  considered  for  brief  analysis  of 
thermal  and  hydraulic  performance  in  the  following  sub-sections. 
Design  and  construction  details  of  such  type  of  conventional 
system  are  described  by  Garg  and  Prakash  [14]. 

2.3.  Thermal  performance 

Thermal  performance  of  a  solar  air  heater  can  be  computed 
with  the  help  of  Hottel-Whillier-Bliss  equation  reported  by 
Duffie  and  Beckman  [15]. 

Qu  =  Ac  Fr[  I(Toc)e—UL(Ti—Ta)  ]  (1) 

or 

Qu  =  =  Fr  [  I(nx)e-UL(Ti-Ta)  ]  (2) 

The  rate  of  valuable  energy  gain  by  flowing  air  in  the  duct  of  a 
solar  air  heater  can  also  be  calculated  from  the  following  equation: 

Q.u=rn  Cp  (To-Tj)  =  hAc  (Tpm-Tam)  (3) 
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Fig.  1.  Conventional  solar  air  heater. 
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The  value  of  heat  transfer  coefficient  ( h )  can  be  increased  by 
various  active  and  passive  augmentation  techniques.  It  can  be 
represented  in  non-dimensional  form  of  Nusselt  number  ( Nu ). 

Nu  =  hl/k  (4) 

Further,  thermal  efficiency  of  a  solar  air  heater  can  be  expressed 
by  the  following  equation; 

>1th  =  y  =  fr  [  (TOC )e-UL(Ti-Ta)/I  ]  (5) 

2.2.  Hydraulic  performance 

Hydraulic  performance  of  a  solar  air  heater  concerns  with 
pressure  drop  (A P)  in  the  duct.  Pressure  drop  accounts  for  energy 
consumption  by  blower  to  propel  air  through  the  duct.  The 
pressure  drop  for  fully  developed  turbulent  flow  through  duct 
with  Re  <  50,000  is  given  as  follows  equation: 


3.2.  Theoretical  approach 

A  theoretical  prediction  works  out  the  consequences  of  a 
mathematical  model,  rather  than  those  of  an  actual  physical  model. 
For  the  physical  processes  of  interest,  the  mathematical  model 
mainly  consists  of  a  set  of  differential  equations.  If  the  methods  of 
classical  mathematics  were  to  be  used  for  solving  these  equations, 
there  would  be  little  hope  of  predicting  many  phenomena  of  prac¬ 
tical  interest.  In  the  theoretical  approach  simplifying  assumptions 
are  used  in  order  to  make  the  problems  tractable  [17], 

3.3.  Computational  approach 

Computational  fluid  dynamics  or  CFD  is  the  analysis  of 
systems  involving  fluid  flow,  heat  transfer  and  associated  phe¬ 
nomena  such  as  chemical  reactions  by  means  of  computer-based 
simulation.  The  technique  is  very  powerful  and  spans  a  wide 
range  of  industrial  and  non-industrial  application  areas.  Some 
examples  are: 


where 

/  =  0.079  RJ°  25  (7) 

2.3.  Thermo-hydraulic  performance 


It  is  necessary  that  while  evaluating  the  performance  of  a  solar 
air  heater  with  respect  to  the  enhancement  of  thermal  gain,  the 
energy  spent  in  propelling  air  should  also  be  taken  into  account. 
It  is  desirable  that  design  of  solar  air  heater  should  be  made 
in  such  a  way  that  it  should  transfer  maximum  heat  energy  to 
the  flowing  fluid  with  minimum  consumption  of  blower  energy. 
Therefore  in  order  to  analyze  overall  performance  of  a  solar  air 
heater,  thermo-hydraulic  performance  should  be  evaluated  by 
considering  thermal  and  hydraulic  characteristics  of  the  collector 
simultaneously.  Thermo-hydraulic  performance  of  a  solar  air 
heater  is  given  by  the  following  index 


Thermohydraulic  performance 


Nu/Nus 

(f/fs)V3 


(8) 


3.  Method  of  prediction 


There  are  basically  three  approaches  or  methods  that  can  be 
used  to  solve  a  problem  of  fluid  flow  and  heat  transfer.  These 
approaches  are: 


1.  Experimental 

2.  Theoretical 

3.  Computational  (CFD) 


3.1.  Experimental  approach 

The  most  reliable  information  about  a  physical  process  is  often 
given  by  actual  measurement.  An  experimental  approach  involving 
full-scale  equipment  can  be  used  to  predict  how  identical  copies 
of  the  equipment  would  perform  under  the  same  conditions.  Such 
full  scale  tests  are,  in  most  cases,  prohibitively  expensive  and  often 
impossible.  The  alternative  then  is  to  perform  experiments  on  small- 
scale  models.  The  resulting  information  however  must  be  extra¬ 
polated  to  full  scale,  and  general  rules  for  doing  this  are  often 
unavailable.  Further,  the  small-scale  models  do  not  always  simulate 
all  the  features  of  the  full-scale  equipment;  frequently,  important 
features  such  as  combustion  or  boiling  are  omitted  from  the  model 
tests.  This  further  reduces  the  usefulness  of  the  test  results.  Finally  it 
must  be  remembered  that  there  are  serious  difficulties  of  measure¬ 
ment  in  many  situations,  and  that  the  measuring  instruments  are  not 
free  from  errors  [16]. 


•  aerodynamics  of  aircraft  and  vehicles:  lift  and  drag 

•  hydrodynamics  of  ships 

•  power  plant:  combustion  in  internal  combustion  engines  and 
gas  turbines 

•  turbo  machinery:  flows  inside  rotating  passages,  diffusers  etc. 

•  electrical  and  electronic  engineering:  cooling  of  equipment 
including  microcircuits 

•  chemical  process  engineering:  mixing  and  separation,  polymer 
moulding 

•  external  and  internal  environment  of  buildings:  wind  loading 
and  heating/ventilation 

•  marine  engineering:  loads  on  off-shore  structures 

•  environmental  engineering:  distribution  of  pollutants  and 
effluents 

•  hydrology  and  oceanography:  flows  in  rivers,  estuaries,  oceans 

•  meteorology:  weather  prediction 

•  biomedical  engineering:  blood  flows  through  arteries  and  veins 

The  equations  governing  the  fluid  flow  are  the  continuity 
(conservation  of  mass),  the  Navier-Stokes  (balance  of  momen¬ 
tum),  and  the  energy  (conservation  of  energy)  equations.  These 
equations  form  a  system  of  coupled  non-linear  partial  differential 
equations  (PDEs).  Because  of  the  coupled  nature  of  the  equations 
and  the  presence  of  non-linear  terms,  the  fluid  flow  equations  are 
generally  not  amenable  to  analytical  method  of  obtaining  the 
solution.  In  general,  closed  form  analytical  solutions  are  possible 
only  if  these  PDEs  can  be  made  linear,  either  because  non-linear 
terms  naturally  drop  out  (as  in  the  case  of  parallel  flows  or  flows 
that  are  inviscid  and  irrotational  everywhere)  or  because  the 
nonlinear  terms  are  small  compared  to  other  terms  so  that  they 
can  be  neglected  (e.g.,  creeping  flows,  small  amplitude  sloshing  of 
liquid  etc.).  If  the  non-linearities  in  the  governing  PDEs  cannot  be 
neglected,  which  is  often  the  case  for  most  engineering  flows,  one 
normally  has  to  resort  to  numerical  methods  to  obtain  solutions. 

CFD  is  the  art  of  replacing  the  differential  equations  governing 
the  Fluid  Flow,  with  a  set  of  algebraic  equations  (this  process  is 
called  the  discretization),  which  in  turn  can  be  solved  with  the  aid 
of  a  digital  computer  to  get  an  approximate  solution.  The  com¬ 
monly  used  discretization  methods  in  CFD  analysis  are  the  finite 
difference  method  (FDM),  the  finite  volume  method  (FVM),  the 
finite  element  method  (FEM),  and  the  boundary  element  method 
(BEM).  Some  special  cases  of  flow  problems  can  also  be  solved 
using  nonstandard  method  like,  boundary  integral  methods,  spec¬ 
tral  methods,  and  pseudo-spectral  methods.  From  the  1960s 
onwards  the  aerospace  industry  has  integrated  CFD  techniques 
into  the  design,  R&D  and  manufacture  of  aircraft  and  jet  engines. 
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More  recently  the  methods  have  been  applied  to  the  design  of 
internal  combustion  engines,  combustion  chambers  of  gas  turbines 
and  furnaces.  Furthermore,  motor  vehicle  manufacturers  now 
routinely  predict  drag  forces,  under-bonnet  air  flows  and  the  in- 
car  environment  with  CFD.  Increasingly  CFD  is  becoming  a  vital 
component  in  the  design  of  industrial  products  and  processes.  The 
ultimate  aim  of  developments  in  the  CFD  field  is  to  provide  a 
capability  comparable  with  other  CAE  (computer-aided  engineer¬ 
ing)  tools  such  as  stress  analysis  codes.  The  main  reason  why  CFD 
has  lagged  behind  is  the  tremendous  complexity  of  the  underlying 
behaviour,  which  precludes  a  description  of  fluid  flows  that  is  at 
the  same  time  economical  and  sufficiently  complete.  The  avail¬ 
ability  of  affordable  high-performance  computing  hardware  and 
the  introduction  of  user-friendly  interfaces  have  led  to  a  recent 
upsurge  of  interest,  and  CFD  has  entered  into  the  wider  industrial 
community  since  the  1990s. 

There  are  several  unique  advantages  of  CFD  over  experiment- 
based  approaches  to  fluid  systems  design: 

•  substantial  reduction  of  lead  times  and  costs  of  new  designs 

•  ability  to  study  systems  where  controlled  experiments  are 
difficult  or  impossible  to  perform  (e.g.,  very  large  systems) 

•  ability  to  study  systems  under  hazardous  conditions  at  and 
beyond  their  normal  performance  limits  (e.g.,  safety  studies 
and  accident  scenarios) 

•  practically  unlimited  level  of  detail  of  results. 

CFD  codes  are  structured  around  the  numerical  algorithms 
that  can  tackle  fluid  flow  problems.  In  order  to  provide  easy 
access  to  their  solving  power  all  commercial  CFD  packages 
include  sophisticated  user  interfaces  to  input  problem  parameters 
and  to  examine  the  results.  Hence  all  codes  contain  three  main 
elements:  (i)  a  pre-processor,  (ii)  a  solver  and  (iii)  a  post¬ 
processor.  We  briefly  examine  the  function  of  each  of  these 
elements  within  the  context  of  a  CFD  code. 

3.3  A.  Pre-processor 

Pre-processing  consists  of  the  input  of  a  flow  problem  to  a  CFD 
program  by  means  of  an  operator-friendly  interface  and  the  subse¬ 
quent  transformation  of  this  input  into  a  form  suitable  for  use  by  the 
solver.  The  user  activities  at  the  pre-processing  stage  involve: 

•  Definition  of  the  geometry  of  the  region  of  interest:  the 
computational  domain 

•  Grid  generation — the  sub-division  of  the  domain  into  a  num¬ 
ber  of  smaller,  non-overlapping  sub-domains:  a  grid  (or  mesh) 
of  cells  (or  control  volumes  or  elements) 

•  Selection  of  the  physical  and  chemical  phenomena  that  need 
to  be  modelled 

•  Definition  of  fluid  properties 

•  Specification  of  appropriate  boundary  conditions  at  cells 
which  coincide  with  or  touch  the  domain  boundary 

Solution  to  a  flow  problem  (velocity,  pressure,  temperature  etc.) 
is  defined  at  nodes  inside  each  cell.  The  accuracy  of  a  CFD  solution 
is  governed  by  the  number  of  cells  in  the  grid.  In  general,  the  larger 
the  number  of  cells,  the  better  the  solution  accuracy.  Both  the 
accuracy  of  a  solution  and  its  cost  in  terms  of  necessary  computer 
hardware  and  calculation  time  are  dependent  on  the  fineness  of 
the  grid.  Optimal  meshes  are  often  non-uniform:  finer  in  areas 
where  large  variations  occur  from  point  to  point  and  coarser  in 
regions  with  relatively  little  change.  Efforts  are  under  way  to 
develop  CFD  codes  with  a  self-adaptive  meshing  capability.  Ulti¬ 
mately  such  programs  will  automatically  refine  the  grid  in  areas  of 
rapid  variations.  A  substantial  amount  of  basic  development  work 


still  needs  to  be  done  before  these  techniques  are  robust  enough  to 
be  incorporated  into  commercial  CFD  codes.  At  present  it  is  still  up 
to  the  skills  of  the  CFD  user  to  design  a  grid  that  is  a  suitable 
compromise  between  desired  accuracy  and  solution  cost. 

33.2.  Solver 

There  are  three  distinct  streams  of  numerical  solution  techni¬ 
ques:  finite  difference,  finite  element  and  spectral  methods.  We  are 
solely  concerned  with  the  finite  volume  method,  a  special  finite 
difference  formulation  that  is  central  to  the  most  well-established 
CFD  codes:  CFX/ANSYS,  FLUENT,  PHOEN1CS  and  STAR-CD.  In  outline 
the  numerical  algorithm  consists  of  the  following  steps: 

•  Integration  of  the  governing  equations  of  fluid  flow  over  all  the 
(finite)  control  volumes  of  the  domain 

•  Discretization — conversion  of  the  resulting  integral  equations 
into  a  system  of  algebraic  equations 

•  Solution  of  the  algebraic  equations  by  an  iterative  method 

The  first  step,  the  control  volume  integration,  distinguishes  the 
finite  volume  method  from  all  other  CFD  techniques.  The  result¬ 
ing  statements  express  the  exact  conservation  of  relevant  proper¬ 
ties  for  each  finite  size  cell.  This  clear  relationship  between  the 
numerical  algorithm  and  the  underlying  physical  conservation 
principle  forms  one  of  the  main  attractions  of  the  finite  volume 
method  and  makes  its  concepts  much  simpler  to  understand  by 
engineers  than  the  finite  element  and  spectral  methods.  The 
conservation  of  a  general  flow  variable  q>,  e.g.,  a  velocity  compo¬ 
nent  or  enthalpy,  within  a  finite  control  volume  can  be  expressed 
as  a  balance  between  the  various  processes  tending  to  increase  or 
decrease  it.  In  words  we  have: 


Rate  of  change 

Net  Rate  of  increase 

of  A>  in  the  control 

of  due  to  convection 

volume  with  respect 

in  to  the  control 

to  time 

volume 

Net  Rate  of  increase  " 

of  <P  due  to  diffusion 

+ 

"Net  Rate  of  creation" 

of  c P  inside  the 

in  to  the  control 

control  volume 

volume 

CFD  codes  contain  discretization  techniques  suitable  for  the 
treatment  of  the  key  transport  phenomena,  convection  (transport 
due  to  fluid  flow)  and  diffusion  (transport  due  to  variations  of  <?>  from 
point  to  point)  as  well  as  for  the  source  terms  (associated  with  the 
creation  or  destruction  of  <?)  and  the  rate  of  change  with  respect  to 
time.  The  underlying  physical  phenomena  are  complex  and  non¬ 
linear  so  an  iterative  solution  approach  is  required.  The  most  popular 
solution  procedures  are  by  the  TDMA  (tri-diagonal  matrix  algorithm) 
line-by-line  solver  of  the  algebraic  equations  and  the  SIMPLE 
algorithm  to  ensure  correct  linkage  between  pressure  and  velocity. 
Commercial  codes  may  also  give  the  user  a  selection  of  further,  more 
recent,  techniques  such  as  Gauss-Seidel  point  iterative  techniques 
with  multigrid  accelerators  and  conjugate  gradient  methods. 

3.3.3.  Post-processor 

As  in  pre-processing,  a  huge  amount  of  development  work  has 
recently  taken  place  in  the  post-processing  field.  Due  to  the 
increased  popularity  of  engineering  workstations,  many  of  which 
have  outstanding  graphics  capabilities,  the  leading  CFD  packages  are 
now  equipped  with  versatile  data  visualization  tools.  These  include: 

•  Domain  geometry  and  grid  display 

•  Vector  plots 

•  Line  and  shaded  contour  plots 

•  2D  and  3D  surface  plots 
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•  Particle  tracking 

•  View  manipulation  (translation,  rotation,  scaling  etc.) 

•  Color  PostScript  output 

More  recently  these  facilities  may  also  include  animation  for 
dynamic  result  display,  and  in  addition  to  graphics  all  codes 
produce  trusty  alphanumeric  output  and  have  data  export  facil¬ 
ities  for  further  manipulation  external  to  the  code.  As  in  many 
other  branches  of  CAE,  the  graphics  output  capabilities  of  CFD 
codes  have  revolutionized  the  communication  of  ideas  to  the  non¬ 
specialist  [18].  Table  1  shows  the  comparison  of  three  approaches. 


4.  Applications  of  CFD  in  various  aspects  of  solar  air  heaters 

A  lot  of  experimental  and  theoretical  studies  have  been 
reported  to  evaluate  performance  of  solar  air  heater.  Kumar 
et  al.  [19]  experimentally  investigated  heat  transfer  and  friction 
characteristics  of  solar  air  heater  by  using  discrete  W-shaped 
roughness  on  one  broad  wall  of  solar  air  heater.  The  maximum 
enhancement  of  Nusselt  number  and  friction  factor  as  a  result  of 
providing  artificial  roughness  was  found  to  be  2.16  and  2.75  times 
that  of  smooth  duct.  Mittal  et  al.  [20]  presented  a  comparison  of 
effective  efficiency  of  solar  air  heaters  having  different  types  of 
geometry  of  roughness  elements  on  the  absorber  plate.  The 
effective  efficiency  was  computed  by  using  the  correlations  for 
heat  transfer  and  friction  factor  developed  by  various  investiga¬ 
tors  within  the  investigated  range  of  operating  and  system 
parameters.  Prasad  and  Saini  [21]  developed  various  design 
curves  for  artificially  roughened  solar  air  heater  that  gave  the 
optimal  thermo-hydraulic  performance.  Prasad  and  Saini  [22] 
investigated  the  effect  of  relative  roughness  pitch  (P/e)  and 
relative  roughness  height  (e/D)  on  the  heat  transfer  coefficient 
and  friction  factor  for  fully  developed  turbulent  flow  in  a  solar  air 
heater  duct  with  small  diameter  protrusion  wires  on  the  absorber 
plate.  It  was  found  that  for  a  given  relative  roughness  pitch,  both 
the  Nusselt  number  and  friction  factor  increased  with  increasing 
relative  roughness  height  and  for  a  given  relative  roughness 
height  both  the  Nusselt  number  and  friction  factor  decreased 
with  increasing  relative  roughness  pitch,  but  not  in  direct  propor¬ 
tion.  Aharwal  et  al.  [23]  carried  out  an  experimental  investigation 
of  heat  transfer  and  friction  factor  characteristics  of  a  rectangular 
duct  roughened  with  repeated  square  cross-section  split-rib  with 


Table  1 

Comparison  of  different  approaches. 


Approach 

Advantage 

Disadvantage 

Experimental 

1.  Capable  of  being  most  realistic 

1.  Equipment  required 

2.  Scaling  problems 

3.  Tunnel  corrections 

4.  Measurement 
difficulties 

5.  Operating  costs 

Theoretical 

1 .  Clean,  general  information, 
which  is  usually  in  formula  form 

1 .  Restricted  to  simple 
geometry  and 
physics 

2.  Usually  restricted  to 
linear  problems 

Computational 

1.  No  restriction  to  linearity 

2.  Complicated  physics  can  be 
treated 

3.  Time  evolution  of  flow  can  be 
obtained 

1.  Truncation  errors 

2.  Boundary  conditions 
problems 

3.  Computer  costs 

a  gap,  on  one  broad  wall  arranged  at  an  inclination  with  respect  to 
the  flow  direction.  The  maximum  enhancement  in  Nusselt  num¬ 
ber  and  friction  factor  was  observed  to  be  2.59  and  2.87  times  of 
that  of  the  smooth  duct  respectively.  The  thermo-hydraulic 
performance  parameter  was  found  to  be  maximum  for  the 
relative  gap  width  of  1.0  and  the  relative  gap  position  of  0.25. 
Multiwork  [24]  carried  out  an  experimental  analysis  of  a  solar  air 
heater  having  V-shaped  staggered  discrete  ribs  on  the  absorber 
plate  and  reported  that  maximum  heat  transfer  enhancement 
occurred  at  an  angle  of  attack  of  60°.  Prasad  and  Mullick  [25] 
utilized  artificial  roughness  in  a  solar  air  heater  duct  in  the  form 
of  small  diameter  wires  to  increase  the  heat  transfer  coefficient 
for  relative  roughness  height  and  relative  roughness  pitch  of 
0.019  and  12.7,  respectively.  Gupta  [26]  investigated  the  effect 
of  relative  roughness  height,  angle  of  attack  and  Reynolds  number 
on  heat  transfer  and  friction  factor  in  rectangular  duct  having 
circular  wire  ribs  on  the  absorber  plate.  The  maximum  enhance¬ 
ment  of  Nusselt  number  and  friction  factor  as  a  result  of  providing 
artificial  roughness  was  found  to  be  1.8  and  2.7  times  that  of 
smooth  duct.  Verma  and  Prasad  [27]  investigated  the  effect  of 
geometrical  parameters  of  circular  wire  ribs  on  heat  transfer  and 
friction  factor.  It  was  observed  that  the  value  of  heat  transfer 
enhancement  factor  varies  from  1.25  to  2.08  within  the  range  of 
parameters  investigated.  The  value  of  optimal  thermo-hydraulic 
performance  was  found  to  be  about  71%  corresponding  to  rough¬ 
ness  Reynolds  number  of  24.  Karwa  [28]  carried  out  a  compara¬ 
tive  experimental  study  of  augmented  heat  transfer  and  friction 
in  a  rectangular  duct  of  a  solar  air  heater  with  rectangular  cross- 
section  ribs  arranged  in  V-continuous  and  V-discrete  pattern. 
The  performance  of  V-down  ribs  was  observed  to  be  better  than 
that  of  V-up  ribs.  Sahu  and  Bhagoria  [29]  experimentally  inves¬ 
tigated  the  effect  of  90°  broken  ribs  as  roughness  elements  and 
found  that  thermal  efficiency  lies  in  the  range  of  51-83.5%.  The 
maximum  enhancement  in  heat  transfer  was  reported  at  the  pitch 
of  20  mm.  Gupta  et  al.  [30]  investigated  the  effect  of  transverse 
and  inclined  wire  roughness  on  fluid  flow  characteristics  of  a 
solar  air  heater.  Authors  reported  that  heat  transfer  coefficient  of 
roughened  duct  using  transverse  wires  fixed  on  absorber  plate  as 
roughness  geometry  can  be  improved  by  a  factor  of  1.8  times  of 
the  smooth  duct  at  the  cost  of  friction  penalty.  Saini  and  Saini  [31  ] 
used  expanded  metal  mesh  as  roughness  geometry  and  obtained 
an  enhancement  of  heat  transfer  coefficient  and  friction  factor  of 
the  order  4  and  5  times  over  the  smooth  duct  corresponding  to  an 
angle  of  attack  of  61.9°  and  72°  respectively.  Momin  et  al.  [32] 
experimentally  investigated  the  effect  of  geometrical  parameters 
of  V-shaped  ribs  on  heat  transfer  and  fluid  flow  characteristics  of 
rectangular  duct  of  solar  air  heater.  The  maximum  enhancement 
of  Nusselt  number  and  friction  factor  as  a  result  of  providing 
artificial  roughness  was  found  to  be  2.3  and  2.83  times  that  of 
smooth  duct.  Karwa  et  al.  [33]  performed  an  experimental 
investigation  of  heat  transfer  and  friction  in  a  high-aspect-ratio 
rectangular  duct  with  repeated  rectangular  cross-section  ribs  on 
one  broad  wall  in  v-discrete  and  discontinuous  patterns  at 
different  angles.  The  performance  of  V-down  ribs  was  observed 
to  be  better  than  that  of  V-up  ribs.  Karwa  et  al.  [34]  carried  out 
an  experimental  investigation  of  heat  transfer  and  friction  of 
rectangular  solar  air  heater  ducts  roughened  with  chamfered 
ribs.  The  highest  heat  transfer  rate  and  also  the  highest  friction 
factor  occurred  for  15°  chamfered  ribs.  The  enhancement  in  the 
Stanton  number  was  reported  to  be  65-90%  while  friction  factor 
was  found  to  be  2.68-2.94  times  as  compared  to  smooth  duct. 
Bhagoria  et  al.  [35]  performed  experiments  to  determine  the 
effect  of  relative  roughness  pitch,  relative  roughness  height  and 
wedge  angle  on  the  heat  transfer  and  friction  factor  in  a  solar  air 
heater  roughened  duct  having  wedge  shaped  rib  roughness. 
Statistical  correlations  for  the  Nusselt  number  and  friction  factor 
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were  developed  in  terms  of  geometrical  parameters  of  the  rough¬ 
ness  elements  and  the  flow  Reynolds  number.  Maximum  heat 
transfer  occurred  for  relative  roughness  pitch  of  7.57.  Wedge 
angle  of  10°  gave  maximum  enhancement  of  heat  transfer.  Saini 
and  Saini  [36]  studied  the  effect  of  arc  shaped  ribs  on  the  heat 
transfer  coefficient  and  friction  factor  of  rectangular  solar  air 
heater  ducts.  Enhancement  of  heat  transfer  and  friction  factor  was 
reported  to  be  of  order  3.6  and  1.75  times  respectively  over 
smooth  solar  air  heater  duct  for  relative  arc  angle  value  of  0.333 
and  relative  roughness  height  value  of  0.042.  Saini  and  Verma 
[37]  investigated  the  effect  of  roughness  and  operating  para¬ 
meters  on  heat  transfer  and  friction  factor  in  a  roughened  duct 
provided  with  dimple-shape  roughness  geometry  for  the  range  of 
Reynolds  number  from  2000  to  12,000,  relative  roughness  height 
from  0.018  to  0.037  and  relative  roughness  pitch  from  8  to  12. 
The  maximum  value  of  the  heat  transfer  was  found  to  be  corre¬ 
sponding  to  relative  roughness  height  of  0.037  and  relative 
roughness  pitch  of  10.  Karmare  and  Tikekar  [38]  performed  an 
experimental  investigation  of  heat  transfer  to  the  airflow  in  the 
rectangular  duct  of  an  aspect  ratio  10:1.  The  top  wall  surface  was 
made  rough  with  metal  ribs  of  circular  cross  section  in  staggered 
manner  to  form  defined  grid.  The  roughened  wall  was  uniformly 
heated  and  the  other  walls  were  insulated.  The  effect  of  relative 
roughness  height  of  grit,  relative  roughness  pitch  of  grit  and 
relative  length  of  grit  on  the  heat  transfer  coefficient  and  friction 
factor  were  investigated  by  authors.  Kumar  et  al.  [39]  carried  out 
an  experimental  investigation  to  determine  the  heat  transfer 
distributions  in  solar  air  heater  having  its  absorber  plate  rough¬ 
ened  with  discrete  w-shaped  ribs.  It  was  reported  that  thermal 
performance  of  the  roughened  duct  was  found  to  be  1.2- 1.8  times 
the  smooth  channel  for  range  of  parameters  investigated.  Jaurker 
et  al.  [40]  investigated  the  effect  of  relative  roughness  pitch, 
relative  roughness  height  and  relative  groove  position  on  the  heat 
transfer  coefficient  and  friction  factor  of  rib-grooved  artificial 
roughness.  The  maximum  heat  transfer  was  obtained  for  a 
relative  roughness  pitch  of  about  6,  and  it  was  decreased  on 
either  side  of  the  relative  roughness  pitch.  The  optimal  condition 
for  heat  transfer  was  found  at  a  groove  position  to  pitch  ratio  of 
0.4  as  compared  to  the  smooth  duct.  As  compared  to  smooth 
surface,  the  presence  of  rib  grooved  artificial  roughness  increases 
the  Nusselt  number  up  to  2.7  times,  while  the  friction  factor 
raises  up  to  3.6  times  in  the  range  of  parameters  investigated. 
Layek  et  al.  [41]  investigated  heat  transfer  and  friction  character¬ 
istics  of  repeated  integral  transverse  chamfered  rib-groove  rough¬ 
ness.  Authors  reported  that  Nusselt  number  and  friction  factor 
was  increased  by  3.24  times  and  3.78  times  respectively  as 
compare  to  smooth  duct.  Maximum  enhancement  of  Nusselt 
number  and  friction  factor  was  obtained  corresponding  to  relative 
groove  position  of  0.4.  Varun  et  al.  [42]  performed  an  experi¬ 
mental  study  on  heat  transfer  and  friction  characteristics  by  using 
a  combination  of  inclined  and  transverse  ribs  on  the  heated  plate 
of  a  solar  air  heater  duct.  The  best  thermal  performance  was 
reported  for  relative  roughness  pitch  value  of  8  and  relative 
roughness  height  value  of  0.030.  Lanjewar  et  al.  [43]  carried  out 
an  experimental  investigation  of  heat  transfer  and  friction  factor 
characteristics  of  rectangular  duct  roughened  with  W-shaped  ribs 
on  its  underside  on  one  broad  wall  arranged  at  an  inclination  with 
respect  to  flow  direction.  Maximum  enhancement  of  Nusselt 
number  and  friction  factor  as  a  result  of  providing  artificial 
roughness  was  found  to  be,  respectively,  2.36  and  2.01  times  that 
of  smooth  duct  for  an  angle  of  attack  of  60°.  Maximum  thermo- 
hydraulic  performance  occurred  at  an  angle  of  attack  of  60°. 
Correlations  were  developed  for  heat  transfer  coefficient  and 
friction  factor  for  roughened  duct.  Varun  et  al.  [44],  Hans  et  al. 
[45],  Bhushan  and  Singh  [46]  and  Kumar  et  al.  [47]  also  reported 
different  investigations  on  roughness  elements  of  different 


shapes,  sizes  and  orientations  carried  out  in  solar  air  heater  in 
order  to  obtain  an  optimum  arrangement  of  roughness  element 
geometry. 

The  literature  survey  reveals  that  a  lot  of  experimental  and 
theoretical  studies  have  been  reported  to  evaluate  performance  of 
solar  air  heater,  but  literature  search  also  reveals  in  this  area  that 
few  CFD  studies  have  been  done  on  the  performance  evaluation  of 
solar  air  heater. 

Arulanandam  et  al.  [48]  investigated  heat  transfer  and  flow  char¬ 
acteristics  for  a  representative  element  of  an  unglazed  transpired- 
plate  absorber  with  holes  in  a  square  pitch  arrangement  under  no 
wind  conditions  using  CFD.  The  study  was  completed  using  TASC 
flow,  a  state-of-the-art  CFD  code.  A  mathematical  model  was  devel¬ 
oped  with  the  relevant  boundary  conditions  and  interfacial  (solid- 
fluid)  conditions  specified.  It  is  reported  that  if  transpired-plate 
absorbers  were  to  be  made  from  lower  conductivity  materials,  such 
as  plastics,  acceptable  efficiencies  could  still  be  achieved.  Fig.  2(a) 
shows  a  sketch  of  the  transpired  plate  with  circular  holes  on  a  square 
pitch  arrangement  in  zero  wind  and  Fig.  2(b)  shows  the  computa¬ 
tional  domain  used  for  CFD  analysis.  The  results  obtained  from  CFD 
approach  were  reported  to  be  closer  to  the  experimental  results. 

Ammari  [49]  proposed  a  CFD  model  for  computing  the  thermal 
performance  of  a  single  pass  flat-plate  solar  air  heater.  Proposed 
solar  air  heater  with  slats  is  shown  schematically  in  Fig.  3. 
Air  channels  were  formed  by  providing  metal  slats  running  along 
the  circulated  air  passage  linking  the  absorber  plate  by  the  bottom 


3  Representative  element 


Fig.  2.  (a)  Sketch  of  the  representative  elements  of  the  absorber  plate  and  (b)  The 
computational  domain. 
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Fig.  4.  (a)  Solution  domain  and  (b)  rectangular  meshing. 


one  in  an  endeavor  to  enhance  the  thermal  efficiency  of  the  solar  air 
collector.  A  computer  code  based  on  FORTRAN  was  developed  to 
solve  for  the  governing  energy  equations  to  estimate  the  mean 
temperatures  of  the  collector.  The  effect  of  volume  airflow  rate, 
collector  length,  and  spacing  between  the  absorber  and  bottom 
plates  on  the  thermal  performance  of  the  present  solar  air  heater 
was  investigated.  Predicted  temperatures  and  efficiencies  of  the  two 
air  heaters  were  presented  in  which  the  proposed  model  heater 
with  the  slats  showed  a  noticeable  improvement  in  performance 
over  that  of  the  most  common  type  conventional  heater.  A  numer¬ 
ical  comparison  of  the  present  design  with  the  most  common  type 
of  solar  air  heaters  was  conducted.  The  CFD  results  showed  good 
agreements  with  experimental  results. 

Chaube  et  al.  [50]  performed  a  computational  analysis  of  heat 
transfer  augmentation  and  flow  characteristics  due  to  artificial 
roughness  in  the  form  of  ribs  on  a  broad  heated  wall  of  a  rectangular 
duct  of  a  solar  air  heater  for  turbulent  flow  (Reynolds  number  range 
of  3000-20,000).  A  2D  analysis  of  heat  transfer  and  fluid  flow 
through  an  artificially  roughened  solar  air  heater  was  carried  out 
using  commercially  available  CFD  software,  FLUENT  6.1.  Shear  Stress 
Transport  k-co  turbulence  model  was  selected  by  comparing  the 
predictions  of  different  turbulence  models  with  experimental  results 
available  in  the  literature.  The  results  obtained  from  CFD  approach 
were  reported  to  be  closer  to  the  experimental  results.  The  solution 
domain  and  rectangular  meshing  used  in  CFD  simulation  are  shown 
in  Fig.  4(a)  and  (b).  Based  on  the  computational  results  following 
conclusions  were  drawn:  (1)  It  was  observed  that  the  2D  analysis 
model  itself  yields  results  which  were  closer  to  the  experimental 
ones  as  compared  to  3D  models  (Fig.  5).  The  3D  model  required 
much  higher  memory  and  computational  time  compared  to  2D 
ones.  Hence,  it  was  sufficient  to  employ  a  simpler  2D  model  which 
was  more  economical  with  the  memory  and  computational  time 
requirements.  (2)  In  the  inter-rib  region,  the  model  predicted  well 
near  the  central  high  heat  transfer  area  but  it  under  predicted 
around  ribs.  (3)  The  peaks  of  local  heat  transfer  were  found  at  the 
reattachment  points.  (4)  The  turbulence  intensity  was  found  max¬ 
imum  at  the  peak  of  the  local  heat  transfer  coefficient  in  the  inter- 
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Fig.  5.  Comparison  between  CFD  predictions  with  experimental  results. 

rib  regions.  (5)  The  peaks  of  the  local  heat  transfer  coefficient  and 
the  reattachment  lengths  decreased  stream  wise  in  the  successive 
inter-rib  regions  up  to  3-4  ribs  and  then  the  distribution  became 
periodic  with  peak  values  and  the  reattachment  lengths  remained 
unchanged.  (6)  The  highest  heat  transfer  was  achieved  with  cham¬ 
fered  ribs  but  the  best  performance  index  was  found  with  rectan¬ 
gular  rib  of  size  3x5  mm. 

Chaube  et  al.  [51]  carried  out  a  2D  computational  analysis  to 
assess  the  comparative  performance  of  the  absorber  plate  of  a 
solar  air  heater  with  different  roughness  elements  using  com¬ 
mercial  software  package  FLUENT  6.1.  The  assessment  was  based 
on  heat  transfer  enhancement  with  minimum  pressure  penalty. 
Ten  different  ribs  shape  viz.  rectangular,  square,  chamfered, 
triangular,  semicircle  etc.  (Fig.  6)  were  investigated  at  the  Rey¬ 
nolds  number  range  from  3000  to  20,000,  in  which  solar  air 
heaters  normally  operate.  The  SST  k-m  turbulence  model  was 
selected  by  comparing  the  predictions  of  different  turbulence 
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Fig.  6.  Different  roughness  elements. 


Fig.  7.  Simulation  of  different  flow  rate. 
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models  with  experimental  results  available  in  the  literature. 
Based  on  the  CFD  results  following  conclusions  were  drawn:  (1) 
The  Shear  Stress  Transport  k-to  turbulence  model  predicted  very 
close  results  to  the  experimental  results,  which  yield  confidence 
in  the  predictions  done  by  CFD  analysis.  (2)  In  transitionally  rough 
flow  a  substantial  enhancement  in  heat  transfer  was  found  along 
with  increase  in  friction  factor,  but  in  fully  rough  flow  a  marginal 
enhancement  in  heat  transfer  over  that  in  the  transitionally  rough 
flow  was  found  with  sharp  increase  (more  than  double  in 
comparison  to  transitionally  rough  flow  regime)  in  friction  factor. 
(3)  There  was  a  great  influence  of  reattachment  of  flow  on 
convective  heat  transfer  coefficient.  (4)  For  rectangular  ribs,  there 
was  an  optimum  width  of  rib  at  which  it  gave  maximum  heat 
transfer  with  minimum  pressure  drop  penalty.  (5)  Chamfered  ribs 
were  not  much  effective  because  although  they  gave  small 
enhancement  in  heat  transfer  but  at  the  cost  of  steep  rise  in 
friction  factor.  (6)  The  thermo-hydraulic  analysis  for  constant 
pumping  power  showed  that  rectangular  rib  (1  mm  x  2  mm)  gave 
the  best  performance  in  comparison  to  other  shapes  under 
consideration. 

Wang  et  al.  [52]  made  comparative  analysis  of  several  types  of 
solar  air  collector  with  numerical  simulation  tools,  got  several 
important  conclusions,  validated  the  advantages  of  unglazed 
transpired  solar  air  collector  UTC  in  ventilation  preheating  field 
and  investigated  the  flow  rate  and  radiant  intensity’s  effort  to  the 
effectiveness.  They  observed  that  CFD  tools  had  its  own  advan¬ 
tages  in  the  field  of  solar  air  collector  research;  it  was  a  powerful 
assistant  tool  for  researching  and  was  worth  to  be  applied  in  more 
abroad  fields.  Based  on  CFD  study  authors  reported:  (1)  Under 
fresh  air  operating  mode,  UTC  had  higher  effectiveness  than  flat 
collector  and  unglazed  collector,  and  could  be  further  improved. 
(2)  As  the  flow  rate  increased,  the  temperature  of  heated  air 
decreased  and  the  effectiveness  increased  (Fig.  7).  (3)  As  the 
radiant  intensity  increased,  the  temperature  of  heated  air 
increased,  but  the  effectiveness  had  almost  no  change. 

Varol  and  Oztop  [53]  carried  out  a  numerical  study  to 
investigate  natural  convection  heat  transfer  and  flow  field  inside 
a  different  shaped  inclined  collectors  such  as  flat  and  wavy 
collector  (Fig.  8).  A  comparison  was  performed  in  these  collectors 
from  the  heat  transfer  point  of  view.  CFDRC  commercial  CFD  code 
was  used  to  simulate  the  laminar  flow  and  thermal  field.  The 
main  purpose  of  this  numerical  work  was  to  compare  the  heat 
transfer  rate  in  flat  and  wavy  collectors.  Based  on  CFD  simulations 
following  conclusions  were  made:  (1)  Natural  convection  heat 


transfer  and  fluid  flow  was  strongly  affected  by  the  shape  and 
inclination  angle  of  the  collector.  (2)  In  all  cases,  more  heat 
transfer  was  obtained  in  the  case  of  wavy  collector  than  that  of 
flat  collector.  (3)  There  could  be  an  optimum  solution  for  the 
governing  parameters  which  was  effective  on  natural  convection 
but  it  was  out  of  the  scope  of  the  study.  (4)  Variation  of  local 
Nusselt  number  was  almost  linear  for  flat  collector.  (5)  With  the 
same  aspect  ratio,  mean  Nusselt  number  increased  with  the 
decreasing  wave  length.  (6)  In  wavy  collectors,  the  mean  Nusselt 
number  decreased  with  the  increasing  inclination  angle.  On  the 
contrary  of  flat  collectors,  the  highest  Rayleigh  number  was 
obtained  at  the  highest  inclination  angle  which  was  the  most 
important  result  of  this  numerical  work. 

Kumar  and  Saini  [54]  carried  out  a  CFD  based  analysis  of  fluid 
flow  and  heat  transfer  characteristics  of  a  solar  air  heaters  having 
roughened  duct  provided  with  artificial  roughness  in  arc  shaped 
geometry  (Fig.  9).  The  heat  transfer  and  flow  analysis  of  the 
chosen  roughness  element  were  carried  out  using  3D  models.  The 
ribs  were  provided  on  the  absorber  plate  whereas  other  sides  of 
the  duct  were  kept  smooth.  FLUENT  6.3.26  commercial  CFD  code 
was  used  for  simulation.  Different  Turbulence  Models  namely 
Shear  Stress  Transport  k-to,  Standard  fc-e,  Renormalization-group 
(RNG)  k-e  and  Realizable  k-e  were  tested  for  smooth  duct  having 
same  cross-section  of  roughened  duct  in  order  to  find  out  the 
validity  of  the  models.  The  results  obtained  by  different  models 
were  compared  with  Dittus-Boelter  empirical  correlation  for 
Nusselt  number.  Based  on  the  computational  results  following 
conclusions  were  drawn:  (1)  Nusselt  number  was  found  to 
increase  with  increase  in  Reynolds  number  whereas  friction 
factor  was  found  to  decrease  with  increase  in  Reynolds  number 
for  all  combinations  of  relative  roughness  height  (e/D)  and 


Fig.  8.  Schematic  models:  (a)  inclined  flat  collector  and  (b)  inclined  wavy  collector. 
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relative  arc  angle  (a/90).  (2)  CFD  results  were  validated  for 
smooth  duct  and  different  CFD  model  results  were  compared 
with  Dittus-Boelter  empirical  relationship  for  smooth  duct.  (3) 
Among  all  the  models  used,  Renormalization-group  (RNG)  k-e 
model  results  were  found  to  have  good  agreement.  (4)  Overall 
enhancement  ratio  with  a  maximum  value  of  1.7  was  found  for 
the  roughness  geometry  corresponding  to  relative  arc  angle 
(a/90)  of  0.333  and  relative  roughness  height  (e/D)  of  0.0426  for 
the  range  of  parameters  considered. 

Karmare  and  Tikekar  [55]  performed  a  CFD  simulation  of  fluid 
flow  and  heat  transfer  in  a  solar  air  heater  duct  with  metal  grit 
ribs  as  roughness  elements  employed  on  one  broad  wail  of  a  solar 
air  heater.  The  circular,  triangular  and  square  shape  rib  grits 
(Fig.  10)  with  the  angle  of  attack  of  54  ;  56°,  58°,  60°  and  62  were 
tested  for  the  same  Reynolds  number.  The  broad  wall  was  subjected 
to  uniform  heat  flux.  The  system  and  operating  parameters  studied 
were:  e/D=0.044,  P/e=  17.5  and  //s=  1.72,  for  the  Reynolds 
number  range  3600-17,000.  To  validate  CFD  results,  experimental 


investigations  were  carried  out  in  the  laboratory.  The  major 
findings  of  the  study  were:  (1)  The  CFD  results  gave  the  good 
agreements  with  experimental  results.  (2)  Artificial  roughness  on 
the  underside  of  the  absorber  plate  of  the  solar  air  heater 
enhanced  the  rate  of  heat  transfer  from  the  absorber  plate  to 
the  flowing  air.  (3)  Amongst  the  different  shape  and  orientations 
analyzed,  the  absorber  plate  of  square  cross-section  rib  with  58° 
angle  of  attack  gave  the  best  possible  results.  (4)  The  percentage 
increase  in  the  heat  transfer  for  58  rib  inclination  plate  over 
smooth  plate  was  found  to  be  about  30%. 

Soi  et  al.  [56]  presented  a  CFD  based  investigation  of  heat 
transfer  and  friction  characteristics  of  artificially  roughened  duct 
of  a  solar  air  heater  (Fig.  11).  K-shaped  (combination  of  transverse 
and  V-up)  roughness  geometry  was  used  as  an  artificial  rough¬ 
ness.  2D  flow  analysis  of  heat  transfer  and  fluid  flow  through  a 
rectangular  duct  roughened  with  K-shaped  ribs  were  carried  out 
using  flow  simulation  module  available  in  SolidWorks  software  of 
version  2010.  2D  computational  domain  and  grid  were  selected. 
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Fig.  10.  Geometry  of  roughened  surface  collector  plate. 
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Fig.  11.  Schematic  of  CFD  model  of  solar  air  heater  duct. 
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Non-uniform  meshing  was  generated  over  the  duct.  FLUENT 
6.3.26  commercial  CFD  code  was  used  for  simulation.  Correlations 
for  Nusselt  number  and  friction  factor  were  also  developed 
corresponding  to  CFD  data  using  SigmaPlot  software.  The  major 
findings  of  the  study  were:  (1)  For  a  given  type  of  artificial 
roughness,  Nusselt  number  increased  with  an  increase  in  Rey¬ 
nolds  number  for  smooth  as  well  as  roughened  plate.  (2)  With  an 
increase  in  relative  roughness  pitch  from  7.5  to  8.3,  Nusselt 
number  increased  and  was  found  maximum  at  P/e  of  8.3.  (3) 
Friction  factor  decreased  with  the  increase  of  Reynolds  number. 
(4)  For  a  given  value  of  relative  roughness  pitch  (P/e)  the  friction 
factor  decreased  with  increase  of  Reynolds  number  5.  A  percen¬ 
tage  deviation  between  CFD  and  predicted  values  of  Nusselt 
number  and  friction  factor  for  K-shaped  ribs  was  found  +15% 
and  +  20%  respectively. 

Giri  [57]  carried  out  a  computational  analysis  of  heat  transfer 
and  flow  characteristics  in  a  solar  air  heater  having  artificial 
roughness  in  the  form  of  different  types  of  ribs  (square,  triangular, 
chamfered  and  semicircular)  on  the  absorber  plate.  CFD  investi¬ 
gations  were  carried  out  for  a  range  of  Reynolds  number  from 
3000  to  15,000  and  relative  roughness  pitch  from  5  to  20.  In  this 
study  CFD  code  FLUENT  6.3.26  was  used  to  perform  a  numerical 
simulation  for  reattachment  point  to  enhance  turbulent  heat 
transfer.  In  this  study  Reynolds-Averaged  Navier-Stokes  analysis 
was  used  as  a  numerical  technique  and  the  RNG  k-e  model  with 
near-wall  treatment  as  a  turbulence  model.  The  results  were 
validated  by  comparing  with  existing  experimental.  Following 
conclusions  were  drawn  from  this  CFD  simulation:  (1)  The  RNG 
k-e  model  gave  good  results  for  the  prediction  of  heat  transfer 
and  friction  characteristics  in  high  aspect  ratio  rib  roughened 
rectangular  duct.  (2)  It  was  observed  that  the  2D  analysis  model 
itself  yields  results,  which  were  closer  to  the  experimental  ones  as 
compared  to  3D  models.  The  3D  model  required  much  higher 
memory  and  computational  time  compared  to  2D  ones.  It  might 
be  because  of  negligible  effect  of  secondary  flow  in  transverse 
ribbed  duct  surface.  Hence,  it  was  sufficient  to  employ  a  simpler 
2D  model  which  was  more  economical  with  the  memory  and 
computational  time  requirement.  (3)  In  the  inter-rib  region,  the 
model  predicted  well  near  the  central  high  heat  transfer  area  but 


it  under  predicted  around  ribs.  (4)  The  peaks  of  local  heat  transfer 
were  found  at  the  reattachment  points.  (5)  The  turbulence 
intensity  was  found  maximum  at  peak  of  the  local  heat  transfer 
coefficient  in  the  inter-rib  regions.  (6)  Chamfered  shape  ribs  gave 
highest  heat  transfer  enhancement.  (7)  CFD  results  also  predicted 
that  Nusselt  number  increased  as  ribs  spacing  increases.  (8) 
Thermal  performance  decreased  as  Reynolds  number  increased. 
(9)  Chamfered  shape  ribs  gave  minimum  friction  as  compared 
to  other  ribs  geometries  (Fig.  12).  (10)  Solar  air  heater  duct 
with  Chamfered  ribs  gave  maximum  heat  transfer  enhancement 
with  minimum  friction  factor  hence  better  thermal  hydraulic 
performance. 

Rajput  [58]  performed  a  CFD  simulation  of  heat  transfer  and 
fluid  flow  behavior  in  a  rectangular  duct  of  solar  air  heater  with 
inverted  U-type  turbulator  roughened  ribs  mounted  on  one  of  the 
principal  wall  (solar  plate).  CFD  code  FLUENT  6.3.26  was  used  to 
perform  a  numerical  simulation.  The  Reynolds-Averaged  Navier- 
Stokes  analysis  was  used  as  a  numerical  technique  and  the  k-e 
turbulence  model  with  near-wall  treatment  as  a  turbulence 
model  (Fig.  13).  The  results  were  validated  by  comparing  with 
existing  experimental  data.  The  major  findings  of  the  study  were: 
(1)  Provision  of  ribs  completely  obstructed  the  viscous  sub-layer 
adjacent  to  the  hot  wall,  which  generated  eddies/recirculation 
zones  upstream  and  downstream,  Eddies  not  only  reduced  the 
heat  transfer  but  also  increased  the  pressure  drop.  With  opening 
the  passage  downstream  the  rib,  intensity  of  eddies  formation  can 
be  reduced.  A  thin  inclined  U-turbulator  geometry  with  a  passage 
underside  can  definitely  prevent  the  formation  of  eddies  and  also 


Fig.  12.  Velocity  profile  for  chamfered  rib. 


Fig.  14.  The  computational  domain. 
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Fig.  13.  Boundary  condition  of  solution  domain. 


72 


AS.  Yadav.J.L.  Bhagoria  /  Renewable  and  Sustainable  Energy  Reviews  23  (2013)  6 0-79 


the  redevelopment  of  two  boundary  layers  at  the  reattachment 
point  in  between  two  adjacent  ribs.  (2)  Roughness  pitch  strongly 
affected  the  flow  pattern  and  also  performance  of  the  duct.  The 
turbulator  geometry  showed  appreciable  heat  transfer  even  at 
low  Reynolds  number  (Re  <  5000)  where  rib  effect  was  inefficient. 
(3)  Value  of  the  Nusselt  number  increased  sharply  at  high 
Reynolds  number  and  this  became  constant  or  increased  very 
slightly  at  low  Reynolds  number.  (4)  In  the  inter-rib  region,  the 
model  predicted  well  near  the  central  high  heat  transfer  area  but 
it  under  predicted  around  ribs.  (5)  Peaks  of  local  heat  transfer 
were  found  at  the  reattachment  points.  Experiments  also  con¬ 
firmed  it.  (6)  Turbulence  intensity  was  found  maximum  at  peak  of 
the  local  heat  transfer  coefficient  in  the  inter-rib  regions.  (7)  The 
turbulences  generated  only  in  the  viscous  sub-layer  region  of 
boundary  layer  resulted  in  better  thermo-hydraulic  performance 
i.e„  maximum  heat  transfer  at  affordable  friction  penalty.  (8) 
Results  obtained  by  Renormalization-group  (RNG)  k-e  model 
were  in  good  agreement  with  the  available  experimental  data. 

Sharma  et  al.  [59]  designed  a  solar  air  heater  duct  in  accordance 
with  the  guidelines  given  in  ASHARE  standard  [1977]  in  order  to 
cany  out  CFD  based  investigation.  CFD  based  performance  analysis 
of  solar  air  heater  duct  provided  with  artificial  roughness  in  the 
form  of  square  type  protrusion  shape  geometry  was  reported. 
2D  computational  domain  and  grid  were  selected.  Non-uniform 
meshing  was  generated  over  the  duct.  Following  conclusions  were 
drawn  from  this  CFD  simulation:  (1)  with  the  application  of  such 
type  of  artificial  roughness,  heat  transfer  coefficient  enhanced  at 
the  cost  of  frictional  penalty.  (2)  Nusselt  number  and  friction  factor 
correlations  were  developed  by  using  CFD  data.  (3)  For  a  given 
value  of  relative  roughness  pitch  (p/e)  Nusselt  number  increased 
with  an  increase  of  Reynolds  number.  Nusselt  number  was  found 
maximum  at  minimum  value  of  relative  roughness  pitch  (p/e).  (4) 
For  a  given  value  of  relative  roughness  pitch  (p/e)  the  friction  factor 
decreased  with  increase  of  Reynolds  number.  (5)  An  Average 
absolute  percentage  deviation  between  CFD  and  predicted  data  of 
Nusselt  number  and  friction  factor  was  found  +  1 5%. 

Sharma  and  Thakur  [60]  conducted  a  CFD  study  to  investigate 
the  heat  transfer  and  friction  loss  characteristics  in  a  solar  air  heater 
having  attachments  of  V-shaped  ribs  roughness  at  60°  relative  to 
flow  direction  pointing  downstream  on  underside  of  the  absorber 
plate.  The  computations  based  on  the  finite  volume  method  with 
the  SIMPLE  algorithm  were  conducted  for  the  air  flow  in  terms  of 
Reynolds  numbers  ranging  from  5000  to  15,000.  Fig.  14  shows  the 
computational  domain  used  for  CFD  analysis.  Figs.  1 5  and  1 6  show 
the  contour  of  heat  flux  on  the  roughened  duct  and  temperature 
contour  of  rib  roughened  plate  in  three  dimensions  obtained  from 
CFD  analysis  respectively.  Following  conclusions  were  drawn  from 
this  CFD  simulation:  (1)  Combined  effect  of  swirling  motion,  detach¬ 
ment  and  reattachment  of  fluid  were  responsible  in  the  increase  of 
heat  transfer  rate  during  CFD  analysis.  (2)  Nusselt  number  increased 
with  increase  in  Reynolds  number  whereas  friction  factor  decreased 
with  increase  in  Reynolds  number  for  all  combinations  of  relative 


Fig.  15.  Contour  plot  of  heat  flux  on  the  roughened  duct. 


Fig.  16.  Temperature  contour  of  rib  roughened  plate. 


roughness  height  (e/D)  and  relative  roughness  pitch  (P/e).  (3)  For 
validating  the  accuracy  of  numerical  solutions,  the  computations 
of  fully  developed  turbulent  flow  forced  convection  in  a  smooth 
rectangular  duct  was  carried  out  to  compare  with  the  exact  solution 
for  Nusselt  number  and  friction  factor,  respectively.  The  present 
numerical  smooth  channel  results  were  found  to  be  in  excellent 
agreement  with  the  exact  solutions  obtained  from  the  literature.  An 
average  absolute  percentage  deviation  between  CFD  and  predicted 
data  of  Nusselt  number  and  friction  factor  was  found  less  than  +  3%. 
(4)  The  v  shaped  rib  roughness  gave  high  rate  of  heat  transfer. 

Gandhi  and  Singh  [61  ]  conducted  a  CFD  study  to  investigate  the 
effect  of  artificial  surface  roughness  on  flow  through  a  rectangular 
duct  having  bottom  wall  roughened  with  repeated  transverse  ribs 
of  wedge  shaped  cross-section.  CFD  investigations  were  carried  out 
for  a  range  of  Reynolds  number  from  5000  to  21,000  for  a  smooth 
duct  and  a  roughened  duct  of  relative  roughness  height  of  0.022, 
relative  roughness  pitch  of  4.5  and  rib  wedge  angle  of  15°. 
Numerical  results  obtained  by  commercial  computational  fluid 
dynamics  (CFD)  code  FLUENT  were  compared  with  the  experi¬ 
mental  results.  Two  dimensional  numerical  modeling  of  the  duct 
flow  using  FLUENT  showed  reasonably  good  agreement  with  the 
experimental  observations  except  for  the  friction  factor.  Following 
conclusions  were  drawn:  (1)  Rib-roughened  surface  in  a  rectan¬ 
gular  duct  produced  much  flatter  velocity  profile  and  higher 
turbulence  intensities  as  compared  to  the  smooth  surface.  This 
enhanced  the  heat  transfer  in  case  of  the  roughened  duct.  (2)  CFD 
analysis  of  the  velocity  profile  and  turbulence  intensity  was  found 
reasonably  well  with  the  experimental  results  for  fc-e  turbulence 
model.  (3)  Friction  factor  for  roughened  duct  was  around  3.4  times 
higher  than  the  smooth  duct  at  Reynolds  number  of  21,000. 

This  review  of  literature  reveals  that  very  few  studies  have 
been  done  on  prediction  of  heat  transfer  and  fluid  flow  in  solar  air 
heater  using  CFD  (Table  2).  CFD  modeling  was  originally  devel¬ 
oped  for  industrial  application.  Today  it  is  being  used  widely  in 
research  work,  product  development,  and  in  almost  all  sphere  of 
activity  where  a  detailed  picture  of  phenomena  involving  heat 
transfer  and  fluid  flow,  phase  change  etc.  is  desired.  Presently  CFD 
techniques  are  increasingly  used  to  model  flow  through  solar  air 
heater.  The  advantages  of  computational  simulations  are  that 
they  can  produce  extremely  large  volumes  of  results  at  virtually 
no  added  expense  and  it  is  very  cheap  to  perform  parametric 
studies  to  optimize  equipment  performance.  The  second  reason 
for  such  work  on  computational  simulation  is  that  some  para¬ 
meters  are  difficult  to  test,  and  experimental  study  is  expensive  as 
well  as  time  consuming.  The  results  obtained  with  the  CFD  are  of 
acceptable  quality. 


5.  Discussion  about  turbulence  models 

The  literature  survey  reveals  that  the  solar  air  heaters  are 
thermo-hydraulically  more  efficient  if  system  operates  at 
Reynolds  numbers  ranges  from  3000  to  19,000.  Reynolds 
number  inside  the  rectangular  duct  of  solar  air  heater  shows 
that  the  flow  is  turbulent.  One  of  the  great  challenges  in  the 
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Table  2 

CFD  simulations  of  solar  air  heater. 


Authors 

Computational  methodology 

Difference  between  experimental  and  simulation  results 

Arulanandam  et  al.  [48] 

•  CFD  code:  TASC  flow 

•  Turbulence  model:  - 

•  Mesh:  uniform 

Close  agreement  observed 

Ammari  [49] 

Chaube  et  al.  [50] 

Fortran 

•  CFD  code:  Fluent  6.1 

•  Turbulence  model:  SST  k-co 

•  Mesh:  rectangular,  non-uniform 

Good  agreement  observed 

Good  agreement  observed 

Chaube  et  al.  [51] 

•  CFD  code:  Fluent  6.1 

•  Turbulence  model:  SST  k-co 

•  Mesh:  rectangular,  non-uniform 

Good  agreement  observed 

Wang  et  al.  [52] 

Varol  and  Oztop  [53] 

•  CFD  code:  CFDRC  ACE+ 

•  Turbulence  model:  - 

•  Mesh:  Non-uniform 

Good  agreement  observed 

Good  agreement  observed 

Kumar  and  Saini  [54] 

•  CFD  Code:  Fluent  6.3.26 

•  Turbulence  model:  RNG  k-s 

•  Mesh:  non-uniform 

Good  agreement  observed 

Karmare  and  Tikekar  [55] 

•  CFD  code:  Fluent  6.2.16 

•  Turbulence  model:  RNG  k-s 

•  Mesh:  non-uniform 

Good  agreement  observed 

Soi  et  al.  [56] 

•  CFD  code:  Fluent 

•  Turbulence  model:  RNG  k-s 

•  Mesh:  non-uniform 

Nusselt  number  ±15% 

Friction  factor  ±  20% 

Giri  [57] 

•  CFD  code:  Fluent 

•  Turbulence  Model:  RNG  k-s 

•  Mesh:  Non-uniform 

Good  agreement  observed 

Rajput  [58] 

•  CFD  code:  Fluent 

•  Turbulent  model:  RNG  k-s 

•  Mesh:  uniform 

Good  agreement  observed 

Sharma  et  al.  [59] 

•  CFD  code:  Fluent 

•  Turbulence  model:  RNG  k-s 

•  Mesh:  non-uniform 

Nusselt  number  ±15% 

Friction  factor  ±15% 

Sharma  and  Thakur  [60] 

•  CFD  code:  Fluent 

•  Turbulent  model:  Realizable  k-s 

•  Mesh:  uniform 

Nusselt  number  ±  3% 

Friction  factor  ±  3% 

Gandhi  and  Singh  [61] 

•  CFD  code:  Fluent 

•  Turbulence  model:  k-co 

•  Mesh:  Uniform 

Good  agreement  observed  except 
for  the  friction  factor 

design  of  a  solar  air  heater  using  CFD  approach  is  the  selection 
of  appropriate  turbulence  model.  When  CFD  simulations  of  a 
flow  through  a  solar  air  heater  are  concerned,  we  often 
wonder,  whether  the  simulation  produces  a  precise  image  of 
the  real  flow.  The  question  is  not  answered  easily,  as  the 
quality  of  the  obtained  result  depends  on  many  parameters  of 
the  computer  model.  Modern  CFD  programs  offer  their  users 
many  possible  parameters  of  the  mathematical  model.  How¬ 
ever,  if  various  turbulence  models  are  available,  the  user  may 
wonder  which  of  them  is  the  best  suited  for  a  solar  air  heater 
system.  The  authors  attempt  to  answer  this  question  in  the 
present  article. 

A  turbulence  model  is  a  computational  procedure  to  close  the 
system  of  mean  flow  equations.  Turbulence  models  allow  the 
calculation  of  the  mean  flow  without  first  calculating  the  full 
time-dependent  flow  field.  Modern  CFD  programs  offer  a  large 


range  of  methods  and  models  to  simulate  turbulence.  The 
methods  include  [62]: 

1.  Reynolds-average  simulation  (RAS)/Reynolds-averaged  Navier- 
Stokes  (RANS):  The  governing  equations  are  solved  in  ensemble- 
averaged  form,  including  appropriate  models  for  the  effect  of 
turbulence.  Different  types  of  incompressible  and  compressible 
RANS  models  are  listed  in  Table  3. 

2.  Large  eddy  simulation  (LES):  Large  turbulent  structures  in  the 
flow  are  resolved  by  the  governing  equations,  while  the  effect 
of  the  sub-grid  scales  (SGS)  are  modeled.  The  scale  separation 
is  obtained  by  applying  a  filter  to  the  governing  equations 
which  also  influences  the  form  of  the  SGS  models.  Different 
types  of  LES  models  are  listed  in  Table  4. 

3.  Detached  eddy  simulation  (DES):  Hybrid  method  that  treats 
near-wall  regions  with  a  RAS  approach  and  the  bulk  flow  with 
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Table  3 

RAS  turbulence  models. 


Table  4 

LES  turbulence  models. 


RAS  turbulence  models 

laminar 

kEpsilon 

kOmega 

kOmegaSST 

RNGkEpsilon 

NonlinearKEShih 

LienCubicKE 

qZeta 

LaunderSharmaKE 

LamBremhorstKE 

LienCubicKELowRe 

LienLeschzinerLowRe 

LRR 

LaunderGibsonRSTM 

realizableKE 

SpalartAllmaras 


for  incompressible  fluids 

Dummy  turbulence  model  for  laminar  flow 

Standard  high-Re  k-e  model 

Standard  high-Re  k-co  model 

k-co-SST  model 

RNG  k-e  model 

Non-linear  Shih  k-e  model 

Lien  cubic  k-e  model 

q-q  model 

Launder-Sharma  low-Re  k-e  model 

Lam-Bremhorst  low-Re  k-e  model 

Lien  cubic  low-Re  k-e  model 

Lien-Leschziner  low-Re  k-e  model 

Launder-Reece-Rodi  RSTM 

Launder-Gibson  RSTM  with  wall-reflection  terms 

Realizable  k-e  model 

Spalart-Allmaras  1-eqn  mixing-length  model 


RAS  turbulence  models  for  compressible  fluids 


laminar 

kEpsilon 

kOmegaSST 

RNGkEpsilon 

LaunderSharmaKE 

LRR 

LaunderGibsonRSTM 

realizableKE 

SpalartAllmaras 


Dummy  turbulence  model  for  laminar  flow 
Standard  k-e  model 
k-co-SST  model 
RNG  k-e  model 

Launder-Sharma  low-Re  k-e  model 
Launder-Reece-Rodi  RSTM 
Launder-Gibson  RSTM 
Realizable  k-e  model 

Spalart-Allmaras  1-eqn  mixing-length  model 


an  LES  approach.  Table  4  also  includes  detached  eddy  simula¬ 
tion  (DES)  models. 

4.  Direct  numerical  simulation  (DNS):  Resolves  all  scales  of 
turbulence  by  solving  the  Navier-Stokes  equations  numeri¬ 
cally  without  any  turbulence  modeling. 

A  number  of  software  based  on  CFD  codes  have  been  devel¬ 
oped,  few  of  them  are:  CFX,  FLUENT,  PHOENICS,  FLOVENT  and 
STAR-CD.  Each  software  is  usually  supported  by  supplementary 
software  for  different  applications,  such  as  domain  model  pre¬ 
paration,  mesh  generation  etc.  Under  the  present  study,  commer¬ 
cial  CFD  code  ANSYS  FLUENT  vl2.1  is  used.  There  are  several 
turbulence  model  included  in  the  commercial  CFD  software 
ANSYS  FLUENT  vl2.1.  The  following  turbulence  models,  available 
in  ANSYS  FLUENT  vl2.1,  are  tested:  Standard  k-&  turbulence 
model,  Realizable  k-s  turbulence  model,  Renormalization-group 
(RNG)  k-s  turbulence  model,  Standard  k-co  turbulence  model, 
and  Shear  Stress  Transport  (SST)  k-co  turbulence  model. 


6.  Selection  of  best  turbulence  model 

In  this  section  a  2-dimentional  CFD  simulation  of  investigated 
in  order  to  find  out  the  best  turbulence  model  for  the  analysis  of  a 
solar  air  heater.  Commercial  CFD  simulation  code  FLUENT  (ver¬ 
sion  12.1)  is  used  as  a  solver.  The  following  assumptions  are 
imposed  for  the  computational  analysis. 

1.  The  flow  is  steady,  fully  developed,  turbulent  and  two 
dimensional. 

2.  The  thermal  conductivity  of  the  duct  wall,  absorber  plate  and 
roughness  material  are  independent  of  temperature. 

3.  The  duct  wall  and  absorber  plate  are  homogeneous  and 
isotropic. 

4.  No-slip  boundary  condition  is  assigned  to  the  walls  in  contact 
with  the  fluid  in  the  model. 

5.  Negligible  radiation  heat  transfer  and  other  heat  losses. 


Isochoric  LES  turbulence  models 

Smagorinsky  Smagorinsky  model 

Smagorinsky2  Smagorinsky  model  with  3D  filter 

dynSmagorinsky  Dynamic  Smagorinsky 

homogenousDynSmagorinsky  Homogeneous  dynamic  Smagorinsky  model 


dynLagrangian 

scaleSimilarity 

mixedSmagorinsky 

dynMixedSmagorinsky 

kOmegaSSTSAS 

oneEqEddy 

dynOneEqEddy 

locDynOneEqEddy 

spectEddyVisc 

LRDDiffStress 

DeardorffDiffStress 

SpalartAllmaras 

SpalartAllmarasDDES 

SpalartAllmarasIDDES 


Lagrangian  two  equation  eddy-viscosity  model 
Scale  similarity  model 
Mixed  Smagorinsky/scale  similarity  model 
Dynamic  mixed  Smagorinsky/scale  similarity 
model 

k-co-SST  scale  adaptive  simulation  (SAS) 
model 

k-e quation  eddy-viscosity  model 
Dynamic  k-equation  eddy-viscosity  model 
Localised  dynamic  k-equation  eddy-viscosity 
model 

Spectral  eddy  viscosity  model 
LRR  differential  stress  model 
Deardorff  differential  stress  model 
Spalart-Allmaras  model 
Spalart-Allmaras  delayed  detached  eddy 
simulation  (DDES)  model 
Spalart-Allmaras  improved  DDES  (IDDES) 
model 


Anisochoric  LES  turbulence  models 

Smagorinsky  Smagorinsky  model 

oneEqEddy  k-equation  eddy-viscosity  model 

dynOneEqEddy  Dynamic  k-equation  eddy-viscosity  model 

lowReOneEqEddy  Low-Re  k-equation  eddy-viscosity  model 

DeardorffDiffStress  Deardorff  differential  stress  model 

SpalartAllmaras  Spalart-Allmaras  1-eqn  mixing-length  model 


LES  deltas 

PrandtlDelta 

cubeRootVolDelta 

maxDeltaxyz 

smoothDelta 


Prandtl  delta 

Cube  root  of  cell  volume  delta 

Maximum  of  x,  y  and  z;  for  structured  hex  cells 

only 

Smoothing  of  delta 


LES  filters 

laplaceFilter 

simpleFilter 

anisotropicFilter 


Laplace  filters 
Simple  filter 
Anisotropic  filter 


The  2D  computational  domain  used  for  CFD  analysis  having  the 
height  of  20  mm,  width  of  200  mm  and  length  of  461  mm.  Fig.  17 
shows  the  2D  computational  domain  used  in  this  study.  Complete 
duct  geometry  is  divided  into  three  sections,  namely,  entrance 
section,  test  section  and  exit  section.  A  short  entrance  length  is 
chosen  because  for  a  roughened  duct,  the  thermally  fully  devel¬ 
oped  flow  is  established  in  a  short  length  2-3  times  hydraulic 
diameter.  The  exit  section  is  used  after  the  test  section  in  order  to 
reduce  the  end  effect  in  the  test  section.  For  the  turbulent  flow 
regime,  ASHRAE  Standard  93-2003  recommends  entrance  and  exit 
length  of  5VWH  and  2.5VWH,  respectively  [63].  The  top  wall 
consists  of  a  0.5  mm  thick  absorber  plate  made  up  of  aluminum. 
A  uniform  heat  flux  of  1000  W/m2  is  considered  for  computational 
analysis.  The  geometrical  and  operating  parameters  employed  in 
this  computational  investigation  are  listed  in  Table  5. 

The  flow  phenomenon  in  a  conventional  solar  air  heater  duct  is 
governed  by  the  continuity,  momentum  and  energy  equations. 
In  the  Cartesian  tensor  system  these  equations  can  be  written  as: 

Continuity  equation: 

i(w)=°  <9) 


Momentum  equation: 


8x, 


,  ,  dP 

«)  =  -Sxi  + 


8xj  f'  \8Xj 


fdUi  i  dUj 

8x  i 


+ 


ox/-^ 


(10) 
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Fig.  17.  Schematic  of  two-dimensional  solution  domain  for  CFD  analysis 


Table  5 

Geometrical  and  operating  parameters  for  computational  analysis. 


Geometrical  and  operating  parameters 

Range 

Entrance  length  of  duct,  Ti’ 

225  mm 

Test  length  of  duct,  ll2' 

115  mm 

Exit  length  of  duct,  ‘Z,3’ 

121  mm 

Width  of  duct,  lW 

100  mm 

Depth  of  duct,  ‘H’ 

20  mm 

Hydraulic  diameter  of  duct,  ‘D’ 

33.33  mm 

Duct  aspect  ratio,  ‘W/H’ 

5 

Uniform  heat  flux,  T 

1,000  W/m2 

Reynolds  number,  ‘Re’ 

3,800-18,000  (6  values) 

Prandtl  number,  ‘Pr’ 

0.7441 

Energy  equation: 


dT 


^(puiri=-((r  +  rt)- 


dXj 


dXj 


(ID 


where  f  and  Ft  are  molecular  thermal  diffusivity  and  turbulent 
thermal  diffusivity,  respectively,  and  are  given  by 


r  =n/Pr  and  rt  =  /j.t/Prt 


(12) 


Table  6 


Thermo-physical  properties  of  working  fluid  (air)  and  absorber  plate  (aluminum) 
for  computational  analysis. 

Properties 

Working  fluid  (air) 

Absorber  plate  (aluminum) 

Density,  ‘ p ’  (kg/m3) 

1.225 

2719 

Specific  heat,  'Cp'  Q/kg//K) 

1006.43 

871 

Viscosity,  (N/m2) 

1.7894e-05 

- 

Thermal  conductivity, 

0.0242 

202.4 

7c’(W/m//K) 


The  following  boundary  conditions  are  assigned  in  FLUENT: 


1.  Velocity  at  inlet  (depending  on  Reynolds  number) 

2.  Turbulence  intensity  at  inlet  (depending  on  Reynolds  number) 

3.  Hydraulic  diameter  of  duct  (fixed) 

4.  Pressure  at  outlet  (fixed  atmospheric  pressure) 

5.  Constant  heat  flux  on  broad  bottom  surface  of  test  section 

6.  No  heat  transfer  from  other  walls  of  the  duct. 

The  properties  of  the  working  fluid  (air)  and  absorber  plate 
material  (aluminum)  are  assumed  to  remain  constant  at  average 
bulk  temperature.  The  thermo-physical  properties  of  working 
fluid  and  absorber  plate  are  illustrated  in  Table  6. 

Meshing  of  the  domain  is  generated  using  ANSYS  ICEM  CFD 
V12.1  software.  A  uniform  mesh  with  very  fine  mesh  size  near  the 
wall  is  used  to  resolve  the  laminar  sub-layer  as  shown  in  Fig.  18. 
Present  uniform  quadrilateral  mesh  contains  45,951  nodes  and 
45,000  cells.  An  exponential  function  is  used  to  concentrate  the 
fine  mesh  near  the  wall  to  reserve  the  high-velocity  gradients 
near  the  walls.  The  first  grid  point  from  the  wall  is  carefully 
adjusted  to  be  located  in  the  linear  region  to  ensure  that  non- 
dimensional  wall  distance  y+  is  less  than  1.0. 

In  the  present  simulation  governing  equations  of  continuity, 
momentum  and  energy  are  solved  by  the  finite  volume  method  in 
the  steady-state  regime  using  the  commercial  CFD  code  ANSYS 
FLUENT  vl2.1.  A  second-order  upwind  scheme  is  chosen  for 
energy  and  momentum  equations.  The  SIMPLE  algorithm  (semi- 
implicit  method  for  pressure  linked  equations)  is  chosen  as 
scheme  to  couple  pressure  and  velocity  [64].  A  uniform  air 
velocity  is  introduced  at  the  inlet  while  a  constant  pressure  outlet 


condition  is  applied  at  the  outlet.  Adiabatic  boundary  condition  is 
implemented  over  the  bottom  duct  wall  while  constant  heat  flux 
condition  is  applied  to  the  upper  duct  wall  of  test  section.  After 
completing  all  the  relevant  settings  FLUENT  starts  to  perform  the 
calculation  in  an  iterative  manner  until  a  sufficient  tolerance, 
defined  by  the  user,  is  achieved.  The  convergence  criteria  of  10_s 
for  the  residuals  of  the  continuity  equation,  10“3  for  the  residuals 
of  the  velocity  components  and  10“6  for  the  residuals  of  the 
energy  are  assumed. 

To  achieve  the  accurate  prediction  of  heat  transfer  and  friction 
factor  in  a  solar  air  heater,  the  predictive  ability  of  five  different 
turbulence  models  including:  The  Standard  k-e  turbulence  model, 
the  Realizable  k-e  turbulence  model,  the  Renormalization-group 
(RNG)  k-e  turbulence  model,  the  Standard  k-co  turbulence  model, 
and  the  Shear  Stress  Transport  (SST)  k-co  turbulence  model,  are 
investigated  and  compared  with  the  available  experimental  data. 

Fig.  19  shows  the  variation  of  average  Nusselt  number  with 
Reynolds  number  for  different  turbulence  models  and  empirical 
correlation  available  for  conventional  solar  air  heater  i.e.,  Dittus- 
Boelter  correlation  [65], 

Dittus-Boelter  equation  :  Nu  =  0.023 Re0SPr04  (13) 

The  average  Nusselt  number  is  observed  to  increase  with 
increase  in  Reynolds  number  for  all  turbulence  models.  It  is 
observed  that  the  results  obtained  by  Renormalization-group 
(RNG)  k-e  model  are  in  good  agreement  with  the  Dittus-Boelter 
empirical  correlation  results.  The  average  absolute  percentage  devi¬ 
ations  between  the  values  predicted  by  Standard  k-co  model  and 
Dittus-Boelter  empirical  correlation  results  is  found  to  be  +  3.58% 
for  Nusselt  number.  Prediction  by  Standard  k-co  model  shows 
more  deviation  with  Dittus-Boelter  empirical  correlation  results. 
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Fig.  19.  Variation  of  Nusselt  number  with  Reynolds  number  using  (a)  Standard  k-e  turbulence  model,  (b)  RNG  k-e  turbulence  model,  (c)  Realizable  k-e  turbulence  model, 
(d)  Standard  k-w  turbulence  model  and  (e)  SST  k-co  turbulence  model. 


The  predictions  with  the  Standard  k-co  model  and  the  SST  k-co 
model  are  over-predicted  in  comparison  with  Dittus-Boelter 
empirical  correlation  results  whereas  Standard  k-e  model  and 
Realizable  k-e,  model  are  under-predicted. 

Fig.  20  shows  the  variation  of  average  friction  factor  with 
Reynolds  number  for  different  turbulence  models  and  empirical 
correlation  available  for  conventional  solar  air  heater  i.e„  Blasius 


correlation  [66], 

Blasius  equation  :/  =  0.0791Re-025  (14) 

It  is  also  observed  that  the  results  obtained  by  Renormalization- 
group  (RNG)  k-e  model  are  in  good  agreement  with  the  Blasius 
equation  empirical  correlation  results.  The  average  absolute  per¬ 
centage  deviations  between  the  values  predicted  by  Standard  k-co 
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Fig.  20.  Variation  of  friction  factor  with  Reynolds  number  using  (a)  Standard  k-e  turbulence  model,  (b)  RNG  k-e.  turbulence  model,  (c)  Realizable  k-e  turbulence  model, 
(d)  Standard  k-w  turbulence  model  and  (e)  SST  k-w  turbulence  model. 


model  and  Blasius  empirical  correlation  results  is  found  to  be 
+  3.28%  for  friction  factor.  Prediction  by  Standard  k-co  model 
shows  more  deviation  with  Blasius  empirical  correlation  results. 

Five  different  turbulence  models,  for  steady  state  conditions, 
available  in  FLUENT,  are  tested  and  the  RNG  k-e  model  is  proven 
as  the  most  appropriate.  The  computed  average  Nusselt  number 
and  friction  factor  are  in  very  good  agreement  with  the  corre¬ 
sponding  experimental  values. 

Finally,  it  can  be  concluded  that  the  Renormalization-group 
(RNG)  k-e  turbulence  model  predicts  very  close  results  to  the 


experimental  results,  which  yields  confidence  in  the  predictions 
done  by  CFD  analysis  in  the  design  of  a  solar  air  heater. 


7.  Conclusions 

This  article  presents  a  detailed  review  of  the  literature  that 
deals  with  the  application  of  CFD  in  design  of  solar  air  heater. 
In  this  article  a  CFD  investigation  is  also  carried  out  to  select  best 
turbulence  model  for  the  design  of  a  solar  air  heater.  A  modern 
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CFD  code  ANSYS  FLUENT  vl2.1  is  used  to  simulate  fluid  flow 
through  a  conventional  solar  air  heater.  A  two-dimensional  flow 
is  assumed.  The  influences  of  the  five  different  turbulence  models 
on  the  quality  of  the  obtained  results  are  tested.  On  the  basis  of 
the  review  of  the  literature  and  CFD  investigation  of  solar  air 
heater,  the  conclusion  can  be  summarized  as  follows: 

1.  This  review  of  literature  reveals  that  a  lot  of  work  has  been 
reported  on  design  of  solar  air  heater  by  experimental 
approach.  This  review  also  reveals  that  a  few  studies  have 
been  done  on  CFD  analysis  of  solar  air  heater. 

2.  CFD  simulation  results  are  found  to  be  in  good  agreement  with 
experimental  results  and  with  the  standard  theoretical 
approaches.  Although  there  are  some  small  discrepancies  due 
to  some  experimental  imperfectness  matters,  we  still  have  a 
good  confidence  in  the  CFD  simulation  program  that  can  be 
used  in  the  future  for  more  complex  solar  air  heater  problem. 

3.  There  is  no  doubt  that  a  major  focus  of  CFD  analysis  of  solar  air 
heater  is  to  enhance  the  design  process  that  deals  with  the 
heat  transfer  and  fluid  flow. 

4.  In  recent  years  CFD  has  been  applied  in  the  design  of  solar  air 
heater.  The  quality  of  the  solutions  obtained  from  CFD  simula¬ 
tions  are  largely  within  the  acceptable  range  proving  that  CFD 
is  an  effective  tool  for  predicting  the  behaviour  and  perfor¬ 
mance  of  a  solar  air  heater. 

5.  One  of  the  great  challenges  in  the  design  of  a  solar  air  heater 
using  CFD  approach  is  the  selection  of  appropriate  turbulence 
model.  The  decision  about  a  suitable  turbulence  model  chosen 
in  a  CFD  computation  is  not  easy. 

6.  The  influences  of  the  five  different  turbulence  models  such  as: 
Standard  k-e  turbulence  model,  the  Realizable  k-e  turbulence 
model,  the  Renormalization-group  (RNG)  k-e  turbulence 
model,  the  Standard  k-co  turbulence  model,  and  the  Shear 
Stress  Transport  (SST)  k-co  turbulence  model  on  the  quality  of 
the  obtained  results  are  tested.  It  appears  from  the  performed 
calculations  that  the  Renormalization-group  k-e  model  yields 
the  best  results  for  two-dimensional  flows  through  a  conven¬ 
tional  solar  air  heater. 

In  summary,  the  purpose  of  this  article  is  to  illustrate  the 
modern  use  of  CFD  in  design  of  solar  air  heater  and  to  use  this  to 
indicate  the  wide-open  future  of  CFD  design.  No  matter  how 
mature  the  techniques  of  CFD  may  become,  the  array  of  future 
and  challenging  applications  of  CFD  is  limitless.  There  is  tremen¬ 
dous  scope  for  future  study  of  solar  air  heater  with  CFD  approach. 
The  information  presented  here  will  be  beneficial  for  beginners  in 
this  area  of  research.  Authors  hope  that  this  article  has  opened  the 
horizons  of  CFD  analysis  of  solar  air  heater  to  researchers. 
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